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Angle-resolved photoelectron spectroscopy is used for a detailed study of the elec- 
tronic structure of the topological insulator Bi2Se3. Nominally stoichiometric and 
calcium-doped samples were investigated. The pristine surface shows the topological 
surface state in the bulk band gap. As time passes, the Dirac point moves to higher 
binding energies, indicating an increasingly strong downward bending of the bands 
near the surface. This time-dependent band bending is related to a contamination 
of the surface and can be accelerated by intentionally exposing the surface to car- 
bon monoxide and other species. For a sufficiently strong band bending, additional 
states appear at the Fermi level. These are interpreted as quantised conduction band 
states. For large band bendings, these states are found to undergo a strong Rashba 
splitting. The formation of quantum well states is also observed for the valence band 
states. Different interpretations of similar data are also discussed. 



INTRODUCTION 



The quantum Hall effect represents a state of matter in which an insulating, two- 
dimensional sample permits perfect conductance along its edges pQ. This is an intriguing 
physical situation with many potential applications, but it can only be realized at low tem- 
peratures and in a high magnetic field. Since its discovery in 1980, researchers have searched 
for this perfect conductance in systems which do not have the constraints of magnetic field 
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and temperature. Twenty five years later, it was suggested that a similar situation could be 
realized using materials with strong spin-orbit splitting of the states without the presence 
of a magnetic field [2H6]. The corresponding state is called the quantum spin Hall effect 
(QSHE) and it was experimentally confirmed in 2007 [7j. In contrast to the quantum Hall 
effect, the QSHE supports counter-propagating states on each sample edge but the propa- 
gation direction is linked to the electron spin. Back-scattering is thus not permitted and, 
again, perfect one-dimensional transport is observed [HOE]. 

While the QSHE corresponds to a two-dimensional system with one-dimensional edge 
states, its three-dimensional analogues are called topological insulators [8HT0]. These mate- 
rials are bulk band insulators but they support metallic surface states. The existence of such 
states is required by the bulk band structure, in contrast to the more coincidental and fragile 
metallic surface states frequently encountered on semiconductor surfaces [llj. The first bulk 
topological insulator found experimentally was obtained by doping the semimetal Bi with 
Sb to obtain Bii_ x Sb x [x ~ 0.1) [T2l - fT4] . a material with a very small band gap, but stable 
surface states that, not surprisingly, strongly resemble those on the pure Bi surface [15j. As 
a material, Bii_^Sb x has several disadvantages for potential applications: its random nature 
(implying poor crystal quality), the very small bulk band gap of less than 30 meV which 
restricts possible applications to low temperatures, and the rather complicated electronic 
structure of its (111) surface. Later, other topological insulators were discovered such as 
Bi 2 Se 3 and Bi 2 Te 3 [T6H20] . These materials are ordered alloys, have a gap size of « 300 meV, 
and their surface electronic structure consists of a single, metallic state with a Dirac cone-like 
dispersion. 

The surface states on a three-dimensional topological insulator show strong similarity 
to the edge states of the QSHE. They are completely non-degenerate with respect to spin, 
except at some symmetry-protected points, and the spin texture of the Dirac cone shows a 
well-defined chirality of the in-plane spin polarisation (201422] . Consequently, the states are 
also immune to back-scattering (23j [24]. Strictly spoken, however, this restriction applies 
to backscattering only, with near back-scattering processes being unlikely but allowed. The 
situation on these surfaces is thus very similar to that for conventional surface states with 
strong Rashba-type spin-orbit splitting (25J [26] and virtually the same as for the surface 
states on the pure semimetal Bi [2TH3T] . 

Recently, the binary alloys Bi 2 Se 3 and Bi 2 Te 3 [19] have been used as prototype topological 
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insulators in experiments because of their simple electronic structure with a single Dirac cone 
at the surface, and because of practical matters such as ease of crystal fabrication and surface 
preparation. The latter is related to the layered structure of Bi 2 Se 3 , as shown in Fig. [jja), 
that is formed of quintuple layers of covalently bonded Bi 2 Se 3 , separated by van der Waals 
gaps. This crystal structure allows samples to be easily cleaved in vacuum thus exposing a 
clean surface. 

The topological character of Bi 2 Se 3 is derived from a parity inversion at the T point 
of the bulk band structure [16]. To see this, consider the so-called time-reversal invariant 
momenta (TRIMs) T[ in the bulk Brillouin zone (BZ), as defined by — I\ = I\ + G, where 
G is a reciprocal lattice vector [32] (see Fig. [jjb)). For each TRIM, it is possible to define 
a parity invariant Si(T{) that accounts for the parity of the occupied bands by taking the 
product of the parity eigenvalues of all (spin degenerate) bands at the TRIM. In Bi 2 Se 3 , 
these parity invariants are all 1, except for the one corresponding to the bulk T point which 
is —1 [16]. Taking the product of the parity invariants at all 8 bulk TRIMs, one thus 
obtains again —1 and this guarantees the bulk material to be a strong topological insulator. 
These parity invariants are then projected onto the surface TRIMs (defined in an analogous 
way), and result in the so-called surface fermion parity. For the hexagonal surface Brillouin 
zone (SBZ) of Bi 2 Se 3 (lll), the T surface TRIM thus obtains a surface fermion parity of 
— 1 while the M TRIM has 1. The surface electronic structure is then closely related to 
these projected parity eigenvalues, with a change in surface fermion parity requiring an odd 
number of Fermi level crossings between the TRIMS and, even more stringently, a surface 
fermion parity eigenvalue of —1 requiring the surface TRIM to be enclosed by an odd number 
of closed Fermi contours [32] . Consistent with this, the T point of Bi 2 Se 3 is enclosed by a 
single, circular Fermi contour. 

Already in the first experimental studies of the topological surface states of Bi 2 Se 3 by 
angle-resolved photoemission (ARPES), a peculiar ageing effect was reported in that the 
observed electronic states gradually move to higher binding energies [dEO]. This effect was 
earlier observed in the similar compound Bi 2 Te 3 [18] and ascribed to band bending induced 
by a slow geometrical relaxation of the van der Waals gap between the quintuple layers of 
the structure. It turns out that this ageing effect brings about a number of interesting and 
controversial spectroscopic observations that will be discussed in this paper. 

The ageing-induced band bending shifts the observed states to higher binding energies 
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FIG. 1: (a) Crystal structure of Bi 2 Se3 and interlayer distances as determined by X-ray diffraction. Three 
primitive lattice vectors are indicated in red while in green the quintuple layer is highlighted, (b) Bulk 
and surface Brillouin zones (BZ and SBZ, respectively) with bulk time-reversal invariant momenta (TRIMs) 
and their projection to surface TRIMs. (c) Bulk band structure along selected high symmetry points and 
projection on the (111) surface after Ref. [33] . 

such that the minimum of the bulk conduction band eventually becomes observable, even for 
initially p-doped samples [HI [20]. Moreover, the spectroscopic appearance of the states near 
the conduction band changes. First, a distinct rim around the bottom of the conduction band 
becomes observable [34], then this rim develops into a well-separated state, and eventually 
it shows a splitting that is strongly reminiscent of a Rashba-type spin-orbit splitting [35] of 
two dimensional states, and is spectroscopically quite similar to Rashba-split surface states 
on surfaces involving heavy elements such as Au or Bi (30], I36H38] . Similar spectroscopic 
features were observed for surfaces which were exposed to magnetic atoms [39j, rest gas in 
the vacuum [40] , carbon monoxide [41] , water [42j and alkali atoms (40], [43J [H] . 

The observed rim around the conduction band was explained in terms of the formation of 
a quantum well state (or two-dimensional electron gas) in the quantum well which is formed 
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because of the strong downward bending of the conduction band at the surface [45]. The 
Rashba-type splitting is a consequence of this quantum well state being placed in a strong 
electric field near the surface [40]. However, alternative explanations have also been given. 
The Rashba-split states have been discussed in terms of being additional, but topologically 
trivial, surface states whose presence becomes necessary because of the strong downward 
band bending [39]. Based on first principles calculations, it has also been suggested that the 
additional states could be caused by the expansion of near-surface van der Waals gaps due 
to the intercalation of adsorbates [461448] . While these calculations show that an increase 
in near-surface van der Waals gaps causes the electronic structure to become more two- 
dimensional, these predictions fail to explain all of the observed changes to the electronic 
structure as will be discussed in more detail below. 

In the present paper we report ARPES results from clean Bi 2 Se 3 surfaces, using crys- 
tals with different bulk doping. We also discuss carbon monoxide adsorption on Ca-doped 
samples. Special emphasis is given to the band-bending induced states in the region of 
the conduction band. The paper is structured as follows: we start with a brief review 
of the ARPES technique and discuss how states of different character (bulk, surface, two- 
dimensional quantum well) are identified in ARPES spectra. We then describe the experi- 
mental and theoretical tools used in this paper. This is followed by a section presenting and 
discussing the results as well as a brief conclusion. 

ANGLE-RESOLVED PHOTOEMISSION 

Even though the ARPES technique is well known and established (for reviews and books 
see Refs. [49H54] ). it is useful to note a few details that will be exploited in the experimental 
investigation presented in the following part. 

We start by briefly reviewing the character of the electronic states that we can expect 
to detect at or near a surface. A schematic summary is given in Fig. [2| Bulk Bloch 
states (Fig. ga)) can always be matched to an exponentially decaying tail outside the 
surface by adjusting the phase between the incoming and reflected wave inside the solid and 
the amplitude of the evanescent wave outside. The bulk electronic structure is thus not 
significantly modified by the introduction of the surface and bulk states. Within (projected) 
band gaps, such as near the BZ boundary in the case of a nearly free electron model, new 
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solutions to the bulk Schrodinger equation can exist. These give real energy eigenvalues, but 
only if a complex k z + iq wave vector is chosen. Clearly, such solutions cannot be normalized 
inside the bulk but at the surface it is in some cases also possible to match them to an 
exponentially decaying wave function in the vacuum. This gives rise to genuine surface states, 
a phenomenon that has been known for a long time [55 , 56j and is textbook material [571158] . 
In the present context, it is worthwhile pointing out that the characteristic exponential 
decay length n into the bulk depends on the energy difference AE between the surface state 
solution and the bulk states, with the order of magnitude given by k ~ ^2mAE/h 2 . For a 
mid-gap state in a topological insulator, this would be quite small, in the order of 1 nm but 
it increases as the state approaches the projected bulk states. Finally, bulk states can be 
localized close to the surface due to quantum-confinement, either by band bending near the 
surface (or interface) of a semiconductor, as routinely realized in semiconductor devices but 
also found at pristine surfaces [59ti63j . or in thin metal films on semiconductors or insulators 
[641 165] . In the case of a semiconductor, the situation is schematically shown in Fig. [2](c). 
The band bending leads to a quantisation of the conduction band states near the surface. 
The wave function for the lowest quantum well state (QWS) is schematically shown. On the 
mesoscopic scale, it is the solution of the Schrodinger equation obtained for a potential well 
of this shape. On a microscopic scale the wave function retains its original Bloch character 
[64H66] . For our purpose, it is worth noticing that the quantum confinement turns the state 
from three-dimensional to two-dimensional in the sense that there is no dispersion in the 
direction perpendicular to the surface. The state is also confined to a region close to the 
surface but its decay away from the surface is not exponential and it typically penetrates 
much deeper into the material than a genuine surface state. Indeed, the physics of these 
quantum-confined states is very different to that of a surface states as they are directly 
derived from bulk states by merely restricting the number of allowed k z values. 

When considering ARPES from the above states, it is useful to first establish the purely 
kinematic conditions for the state's observation. The introduction of the surface breaks 
the translational periodicity of the crystal in the z direction and the wave vector in that 
direction k z is thus no longer well-defined. The components of the wave vector parallel to 
the surface (fcy), on the other hand, are still well defined and must be conserved in the 
photoemission process. They are thus obtained from the photoemission angle with respect 
to the normal of the sample and the kinetic energy of the photoelectron. In order to map 
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FIG. 2: Electronic states near the surface of a solid, (a) Dispersion of nearly free-electron like bulk 
states and the wave function of such a state near the surface (real part). The one-dimensional lattice 
constant perpendicular to the surface is a. The Bloch wave inside the crystal can always be matched to an 
exponentially decaying tail outside the crystal, (b) (Projected) bulk band gaps can give rise to new Bloch 
states with a complex wave number k z + iq, i.e. an exponential growth in one direction. In some cases, these 
exponentially increasing solutions can also be matched to an exponential tail outside the surface, giving rise 
to a surface state, (c) Left: In a semiconductor, a strong downward bending of the conduction band can lead 
to the quantisation of the conduction band states near the surface. The envelope wave function of the lowest 
state (solid line) has no node and the microscopic state still has Bloch- wave character (dashed line). Right: 
In a band structure picture, this quantisation reduces the permitted /^-values in a band (only considered 
for the lower band). The states thus loose their three dimensional dispersion and become two-dimensional. 
The near-surface Bloch wave looks similar to that of the parent bulk state. 
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the photoemission intensity on the hemisphere over the sample's surface, it is possible to 
implement an azimuthal scan (rotating the sample around the surface normal direction) 
for every polar angle considered. In today's ARPES setups, the electron analyzer often 
contains a two-dimensional detector that can simultaneously image the energy dispersion of 
the electrons and the angular dispersion in one direction (see Fig. [3]). In this case, it is useful 
to consider the dispersion along two polar angles (0 and 0), along orthogonal directions as 
presented in Fig. [31(d). In this configuration, the relation between the angles and kn is 



k|| = (sin (0)x + cos (0) sin (0)y) 



2m e E kin 



h 2 



where E^ [n is the kinetic energy of the photoemitted electron. 
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FIG. 3: Scheme depicting the hemispherical analyser and the geometry of the ARPES experiment, (a) 
Cross section and definition of the dispersive plane of the analyser with the orbits of the higher (lower) 
kinetic energy electrons in red (blue), (b) Cross section of the non dispersive plane of the analyser with the 
orbits of the electrons emitted at an angle +/3 (— j3) in green (blue), (c) CCD picture of the acquired spectra 
on the detector, d) Scheme presenting the and <\> angles used to map the photoemission intensity in the 
experiments performed. 



As pointed out above, the wave vector in the z direction is not a good quantum number 
near the surface because the translational symmetry in this direction is broken. Even if we 
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ignore this fundamental problem for the time being, it would not be possible to recover k z 
inside the sample from the value measured outside due to the potential change at the surface 
that leads to a refraction effect at the surface barrier [53]. Despite these difficulties, it is 
often still possible to recover k z of the initial state from the photoemission spectra, provided 
that the dispersion of the final states is known. This is illustrated in Fig. |4](a). As the 
photon energy is changed, different initial state energies and k z values can be probed and 
the bulk state in question will appear as a peak at different binding energies in the spectrum, 
depending on the photon energy used. The state's binding energy at high symmetry points 
( k z = nn j a where n = 0, 1, 2... and a is the one-dimensional lattice constant perpendicular 
to the surface) can easily be read from the spectra, but the detailed dispersion of the state 
can only be recovered if the final state dispersion is known. 

Even if this is not the case, it is often possible to obtain a good approximation of the 
entire dispersion by assuming final states of a simple form, for example free electron final 
states. In this case it is possible to calculate k z plus or minus a reciprocal lattice vector as 
a function of the emission angle 6 and the inner potential Vo using 

K = V 2me /^(^o + ^kincos (9)). (2) 

The inner potential Vo can be determined iteratively by requiring the resulting k z to be 
consistent with the binding energy extrema at symmetry points. 

The problem that k z is not well-defined anymore cannot be neglected. In particular, the 
finite escape depth of the photoelectron leads to a k z broadening 5k z . As the state disperses 
in the k z direction, this also leads to an energy broadening 5E, as illustrated in Fig. |i](a). 

This problem is absent for photoemission from surface states and quantum-confined 
states, shown in Fig. [4](b) and (c), respectively (again, the quantum confinement has only 
been introduced for the lowest band). Strictly spoken, neither the surface state nor the quan- 
tum confined state have a well-defined k z because they are localized in the z direction. This 
is usually a much more severe restriction for the more strongly confined surface states than 
for the deeply penetrating QWS. In neither case does the wave function loose its periodicity 
perpendicular to the surface completely and the Bloch wave character in the z direction is 
partly retained. Thus, the markers in the figure symbolize the k z value the state is derived 
from and the horizontal lines symbolize that k z is not defined anymore. In any case, the 
absence of a dispersion with k z implies directly that a k z broadening does not result in an 
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energy broadening in the photoemission spectra, unlike the case for the bulk states. This 
implies that it is often easier to determine a bulk band structure in the presence of quantum 
confinement than in its absence [64] . The fact that k z is ill-defined also implies that surface 
states and QWS can be observed not only for one photon energy but for a broad range of 
photon energies. 
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FIG. 4: Schematic picture of the photoemission process from different types of electronic states. E F and 
Ey are the Fermi and vacuum level, respectively, (a) Bulk states measured with different photon energies 
appear at different binding energies in the spectrum, i.e. they show dispersion. The fact that k z is not 
well-defined leads to a broadening Sk z that is then reflected in an energy broadening SE of the observed 
peaks, (b) Quantum well states are derived from restricting the possible k z values in a structure (green 
markers, only shown for the lowest band) but the k z broadening (green lines) implies that they can be 
observed for a range of photon energies. As they do not show dispersion, this will not lead to an energy 
broadening. Similar arguments hold for surface states (c). 



With this we come to the dynamic properties of the photoemission from bulk states, 
surface states and QWS, and we consider the photoemission intensity we can expect to 
measure. Even if we assume perfect angular and energy resolution in the experiment, the 
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complete expression is rather complicated [521 • O ne obtains: 

/(£«„, k) oc ^^(k^k^p/^-^fcm-^T) y ^(^-^-^k)^,^)^, (3) 

where and kj are the three-dimensional wave vectors for the initial and final state in the 
solid, k is the wave vector for the free electron outside the solid, hv is the photon energy, 
/ is the Fermi distribution, $ is the work function, A is the hole spectral function and C 
is a Lorentzian distribution to account for the real-space damping of the outgoing electron 
wave, with k® being the perpendicular wave vector for the un-damped final state wave. Mji 
is the matrix element for the photoemission process. The integration is over all possible 
perpendicular wave vector components k z . For QWS or surface states, this integration 
becomes irrelevant. 

The quantity of most interest is often the spectral function A. In simple terms, this is a 
measure of finding an electron with a certain energy and crystal momentum, i.e. an image 
of the band structure, including lifetime broadening effects. For two-dimensional systems 
such as surface states and QWS, the photoemission intensity is directly proportional to A, 
as no final state broadening has to be considered. For a detailed further analysis it is only 
necessary to assume that the matrix element does not depend strongly on the binding energy 
or wave vector. 

In the present context, the matrix element depends rather strongly on the photon energy 
and this dependence contains valuable information on the nature of the observed states. 
This is evident when considering Figs. [2] and [4j For bulk states, no significant change in 
the photoemission intensity is expected as the photon energy is varied (assuming a gap-free 
final state) . A different photon energy will merely lead to emission from a bulk state around 
a different k z (Fig. [4](a)). For surface states and QWS (Fig. [i] (b) and (c)) this is not so: 
both lack a well-defined k z but the wave functions still have a remaining Bloch character 
perpendicular to the surface. For the QWS, this Bloch character is very similar to the 
parent bulk states and for the surface state it is similar to the bands near the BZ boundary 
at k z = n/a (in this particular case). Thus, a resonant enhancement in the photoemission 
intensity of such states is expected for photon energies that would have lead to emission 
from the parent bulk states. This is indeed observed in many cases J67J [68] and it will be 
used in the following discussion. 
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METHODS 

Samples grown from stoichiometric mixutures are known to be highly n-doped due to 
presence of defects in the bulk and mainly charged Se vacancies [69]. For this reason two 
different Bi 2 Se 3 crystals were used during the experiments. The first batch of crystals, 
named as intrinsic, was obtained from a Se-rich mixture of 5N purity elements (Bi : Se = 2 
: 3.3) melted at 860°C for 24 hours in an evacuated quartz ampoule, cooled down to 650°C 
at a rate of 2.5°C/h, and then annealed at 650°C for 7 days. A second batch, named as 
Ca-doped, was obtained by a mixture of Bi (5N purity) , Se (5N purity) and Ca (99.5% 
purity) with ratio Bi : Se : Ca = 1.996 : 3 : 0.004 melted at 860°C for 24 hours in an 
evacuated quartz ampoule. Once cooled down from 860°C to 750°C at a rate of 50°C/h and 
then from 750°C to 600°C at a rate of 2 °C/h, it was annealed at 600°C for 7 days. The 
Ca doping is able to recover the insulating state and the long cooling time produces large 
crystals. The crystal structure and quality was determined using x-ray diffraction at room 
temperature. 

ARPES spectra were acquired at the SGM III undulator beamline at the synchrotron 
light source ASTRID [70J which is equipped with a SPECS Phoibos 150 spectrometer. The 
combined energy resolution (photons and electrons) in this work was better than 15 meV and 
the angular resolution was better than 0.13°. After cleaving the sample in situ and at room 
temperature (pressure better than 1 x 10 -8 mbar), it was immediately transferred to the 
measurement position at a pressure better than 4 x 10 -10 mbar and cooled to a temperature 
of 60 K. 

Simple model calculations were used to determine the energy positions of the QWS in the 
near-surface potential well that is caused by the downward bending of the conduction band. 
To this end, the potential near the surface was simulated using a Schottky model. A constant 
charge density p was assumed to be present between the surface [z = 0) and a certain depth 
{z — A). From this, the electrostatic potential was calculated via the Poisson equation. We 
then numerically solved the Schrodinger equation in this potential and obtain the allowed 
eigentstates and their energies. The model has two parameters, p and the width of the 
charge layer A. These can be determined by two conditions. The first is that the potential 
must reproduce the experimentally observed band bending. The experimental value can be 
inferred from the position of any sharp state in the spectrum. Either a valence state or a 
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core level can be used [71J. We have used the topological surface state for this purpose (see 
below) . The second condition is that the solution of the Schrodinger equation must reproduce 
the number and approximate energies of the experimentally observed occupied QWS. The 
solution of the Schrodinger equation also yields the envelope wave functions as shown in Fig. 
[2|c). As usual with this type of calculation, the Bloch character of the states is taken into 
account by using the correct effective mass of the electrons in the z direction. Here this is 
set to 0.24 m [34]. Note that this simple Schottky model is quite different from a previously 
used model [341 [61] that takes nonparabolic bands into account via the k • p approximation 
and calculates the quantum well states via a coupled solution of the Schrodinger and Poisson 
equations [6]] . The purpose of the calculations here is not a quantitative agreement between 
experiment and model or an extraction of parameters like the near-surface charge density. It 
is merely to have a fast calculation that can produce quantum well states for different degrees 
of band bending and this can be handled satisfactory by this and even simpler models, such 
as a triangular quantum well [42] . 



RESULTS AND DISCUSSION 



We will start by discussing the pristine surface of intrinsic and Ca-doped Bi 2 Se 3 crystals. 
Fig. [5](a) and (b) show ARPES spectra taken from the two batches and in the simple picture 
described above, these can be viewed as cuts through the spectral function along a certain 
direction in the SBZ. Fig. [5](a) shows the the dispersion of the states, i.e. the photoemission 
intensity as a function of k in the KTK direction and binding energy, whereas Fig. [5](b) 
shows the photoemission intensity at the Fermi level as a function of ky. Both spectra are 
cuts through a three-dimensional data set in which the photoemission intensity is collected 
as a function of binding energy and the two components of ky . The distinct V-shaped feature 
in Fig. [5](a) stems from the topological surface state, forming a conical dispersion with its 
Dirac point at a binding energy of 360 meV. At the Fermi level, i.e. far away from the 
Dirac point, the dispersion deviates from a simple cone and the constant energy surfaces 
are hexagonal (Fig. [5] (b)) [34j [72]. The anisotropy of the Fermi velocity is known and 
explained for the analogous Bi 2 Te 3 case where the deviation from a perfect Dirac cone is 
more evident [T9l 173] . The features at the highest binding energies and the region of high 
intensity inside the Dirac cone are due to emission from the bulk valence and conduction 
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bands, respectively. Despite the Se enriched mixture used in the growth, these samples are 
still degenerately n-doped. We also see that the spectral features linked to the topological 
state, a genuine surface state, are considerably sharper than those of the other states. This 
is due to the fc z -smearing explained in Fig. [4j 




FIG. 5: ARPES spectra for the pristine surface of Bi2Se3. High photoemission intensity is displayed in 
bright, (a) Energy dispersion in the KTK direction of the SBZ and (b) Fermi surface for the stoichiometric 
Bi2Se3 sample, (c) and (d) Energy dispersion and Fermi surface for the Ca-doped sample, respectively. The 
Dirac point is at a binding energy of ~ 50 meV. 

The Ca doping is able to counteract the natural doping, resulting in a shift of the Dirac 
point from a binding energy of 365 meV to 47 meV as shown in Fig. [5] (c) and (d). In this 
case, the topological state appears broader, consistent with the higher degree of disorder 
due to the random distribution of dopants in substitutional Bi sites [74] . Such disorder is 
relevant even for surface states, as it induces an additional smearing in k\\. 

The dispersion of the bulk states in the direction perpendicular to the surface can be 
inferred from data taken as a function of the photon energy, such that different values of k z 
are probed (see Fig. [i](a)). Fig. [6] shows the result of such a photon energy scan, presenting 
the photoemission intensity at normal emission acquired for the intrinsic sample. Data were 
taken as a function of photon energy, but the horizontal axis has already been transformed 
to k z values instead of photon energies, using ([2]). For this transformation, an inner potential 
of V =11.8 eV was assumed [171153] . 

The dispersing features stem from the conduction band (CB, blue) and valence band 
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(VB, magenta). The bottom of the conduction band appears together with the top of the 
VB at k z « 2.63 A -1 (hz/ « 19.2 eV), corresponding to the bulk V point. The bottom of 
the VB is reached at k z « 2.96 A -1 (hz/ « 26.6 eV) and corresponds to the bulk Z point. 



The green line in the scan corresponds to the Dirac point of the surface state. The position 
of the Dirac point cannot be easily read from this scan but it is inferred from the dispersion 
of the topological state in images such as Fig. [5] (a). In fact, the data set of Fig. [6] was 
constructed by taking the centre intensity of many of such images taken at different photon 
energies. Note that the position of the Dirac point does not appear at a strictly constant 
energy, even though this should be expected for a two-dimensional state. The reason for 
this is the ageing effect with the time-dependent downward band bending. Taking a data 
set large enough to construct Fig. [6] takes approximately 30 min, and this is long enough 
for the ageing effect to be noticeable. 



FIG. 6: Photon energy scan on the pristine surface of intrinsic Bi2Se3 crystal illustrating the dispersion 
of the states at normal emission as a function of k z . The data shown are a subset of a larger photon energy 
scan between hv — 14 eV and hv = 32 eV. The drift of the Dirac point with photon energy is due to the 
ageing effect that occurs during the scan. The CB and VB (highlighted with blue and magenta lines as a 
guide to the eye) disperse, revealing the bulk T and Z points. The dashed magenta line is a shifted replica 
of the valence band dispersion caused by a surface umklapp process. 



This corresponds to a TZ distance of 0.329 A 1 and is consistent with the literature 
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In addition to the expected valence band dispersion marked by the magenta line, a second 
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identical but weaker dispersion is found (dashed magenta line). This is shifted by half a 
reciprocal lattice vector in the k z direction, i.e. it has its highest binding energy at T instead 
of Z. The presence of this band is attributed to an umklapp process involving a surface 
reciprocal lattice vector. Such processes are well known for similar BZ shapes, where the k z 
value of the T point in the first Brillouin zone is the same as for the zone boundary in the 
neighbouring zone [68J. 

Fig. [7] provides a more detailed picture of the ageing effect. As soon as one hour after 
cleaving the sample, a clear shift of the Dirac point is observable together with the appear- 
ance of additional states: a parabolic state at the rim of the CB and a sharp M-shaped state 
in the VB region. As time elapses the band bending increases and the number of additional 
states increases. Eventually, two distinct states are found the CB region, with the lower 
state showing a clear Rashba splitting [40]. The expected spin-texture of these states has 
been confirmed by spin-resolved ARPES [40j . The ageing also leads to a small broadening of 
the topological state, indicating surface disorder [76]. It should be stressed that the drift of 
the photoemission features towards higher binding energies is a well known phenomenon for 
semiconductors, indicating a downward shift of the electronic bands near the surface [77H80] 

The states in the conduction band have been interpreted as QWS caused by the strong 
band bending near the surface. The situation for the lower conduction band QWS cor- 
responds to that in Fig. [2](c) and its energy position can indeed be reconciled with the 
observed amount of band bending, using either a simple Schottky model (see below) or a 
more sophisticated solution of the coupled Poisson and Schrodinger equations [3H EI] . The 
two-dimensional character of the states can be tested directly by a photon energy scan. 
The result is given in Fig. [7|(f) and should be compared to Fig. |6j The difference in the 
conduction band region is clear: in the case of weak band bending, the dispersion of the en- 
tire conduction band is clearly seen, but for stronger band bending only two non-dispersive 
features are observed at similar energies. This is a clear experimental indication for the for- 
mation of QWS. It is also seen that the intensity of the non-dispersive feature is resonantly 
enhanced at the V point, i.e. at the energy where previously the non-confined state was 
observed. 

In the following we will argue that the band bending is induced by the adsorption of 
impurities. Indications of this are the stronger tendency towards ageing at lower tempera- 
tures and higher background pressures and that the band bending can be partly reversed 
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(a) 



10 min from cleave 



(b) 



1 h from cleave 



(c) 



3h from cleave 



(d) 



16h from cleave 




FIG. 7: ARPES spectra acquired at (a) 10 min, (b) 1 h, (c) 3 h and (d) 16 h after cleave showing 
photoemission intensity as a function of binding energy and fey (hv — 16 eV). Additional states in the CB 
are clearly visible. The band bending is indicated by the shift of the Dirac points towards higher binding 
energy, (e) Magnification of the Rashba split quantum well state of figure (d). (f) Photoemission intensity 
at normal emission as a function of binding energy and k z acquired with a photon energy scan. The data 
shown in (f) are a subset of a larger energy scan between hv = 14 eV and hv = 32 eV. 



by annealing the surface to higher temperatures [40]. We will demonstrate how the same 
ageing effect can be obtained in a controllable and tunable way by depositing known species 
on the surface. 

To test which species can cause the band bending, the clean surfaces have been exposed 
to a controlled dose of carbon monoxide [41], as CO represents one of the common rest gas 
molecules present in an ultra-high vacuum recipient. Similar tests have been performed with 
CO2 and N 2 but this did not reveal an appreciable acceleration of the normal ageing effect 
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(data not shown), apart from an increased broadening and a signal degradation. 

Fig. [8] is analogous to the time dependence shown in Fig. [7| but here spectra were 
acquired at a partial pressure of 8 x 10 -9 mbar of CO on a Ca-doped sample that was kept 
at a temperature of 65 K. Acquisition was started 1 hour after the cleave. The development 
of the band bending is faster and more pronounced than before. Fig. [8j[a)-(c) show the 
dispersion at different times. These spectra are merely a small subset of a large series 
of images. A cut through the centre of all images (i.e. the energy and time-dependent 
photoemission intensity at normal emission) is shown in Fig. [8](d). This illustrates the 
position of the states as the band bending increases. In contrast to Fig. [7[ states of the 
non-confined conduction band are never observed. This is ascribed to the fact that we start 
out with a Ca-doped sample and the conduction band only becomes visible when the band 
bending is already considerable. For a sufficiently large dose of CO, after circa 20 minutes 
of exposure, the band bending is strong enough for the lowest CB QWS to move below 
the Fermi level and thus become observable. In addition to the time of exposure, the band 
bending values are given on the horizontal axis. The values have been extracted from the 
position of the topological state in images like Fig. [8] for a non-zero value of k X) where the 
state is always clearly observable, even for strong band bending. Saturation is reached when 
the band bending is stabilized and the the binding energies of the QWS have reached a 
steady value. The saturation coverage of CO has not been determined but it is expected to 
be < 1 monolayer. Multilayer formation has not been observed and it is not expected on 
inert surfaces at these temperatures and pressures. 

Fig. |8j[d) also shows the results of the Schottky model simulation. The energies of the 
calculated QWS are plotted as blue, azure and cyan dashed lines on top of the data. Clearly, 
the agreement between the simple model and the experiment is very good. As discussed 
previously, the model has two free parameters, the charge density p and the depth below the 
surface A. These parameters can be adjusted to achieve a good fit of the QWS energies but 
here they are only adjusted for one value of the coverage and then kept fixed for all others. 

Again, a photon energy scan confirms the two-dimensional character of the of the con- 
duction band QWS (Fig. [9]). It also shows that the VB states are now quantised, too, with 
the previously continuous valence band now being replaced by distinct features. These are 
also visible Fig. |8j and have a characteristic M-like shape. This simultaneous quantisation 
of both valence band and conduction band states is rather unusual and related to the large 
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(a) (b) (c) 




exposure time (min) 

FIG. 8: ARPES spectra acquired during the controlled dose of CO on Ca-doped Bi2Se3 (hv = 16 eV). The 
acquisition and dosing was started about 1 h after the cleave, (a)-(c) Photoemission intensity as a function 
of binding energy and fey at different exposure times, (d) Photoemission intensity at normal emission as a 
function of binding energy and exposure time. The states in the marked regions are the quantum well states 
(CB QWS and VB QWS). The Dirac point is indicated by the green arrow while the blue, azure and cyan 
dashed lines show the energies of the calculated QWS. 

projected band gap under the top part of the VB shown in Fig. [jjc). It is explained in Ref. 
[4T] and shall not be discussed in much further detail here. However, it is interesting to 
note the enhancement of the M-shaped states for particular photon energies. This happens 
at the same energies as the emission from the corresponding non-quantised valence band 
states through direct photoemission ( in the region of hv « 25.3 eV) and via the umklapp 
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process with a surface reciprocal lattice vector [hv « 19 eV). This behaviour is consistent 
with what would be expected for quantum confined states in the VB and with the simple 
picture given in Fig. [4j 




FIG. 9: Photon energy scan after exposing the Ca doped samples to carbon monoxide, (a) Photoemission 
intensity at normal emission in function of binding energy and k z . (b)-(d) Photoemission intensity at selected 
k z values depicting the energy dispersion relation 

Interestingly, in both this work, and for Bi 2 Te 3 [81j, the QWS are not found to be 
Rashba-split, in contrast to the QWS in the conduction band. The reason for this is likely 
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FIG. 10: Dispersion and Fermi surface map after CO dosing (hv = 16 eV) on Ca doped samples, (a) and (c) 
Photoemission intensity in function of binding energy and fey in the KTK and Mf M directions, respectively, 
(b) Fermi surface map around f showing the hexagonal warping of the QWS that diminishes as the distance 
from f decreases, (d) Spectrum taken at 250 K after CO dosing showing the partial reversibility of the band 
bending. 

the different microscopic wave functions of the VB and CB states. It is well-known that the 
Rashba splitting of extended states is largely an atomic effect and a detailed understanding 
of the size [82] and even the sign [83] of the splitting requires a knowledge of the atomic wave 
function near the nuclei. Indeed, the different wave function character is already evident 
from the state's dispersion near f where the VB states are much flatter and a small Rashba 
splitting, being an offset in fc, would be hard to detect. 



The Fermi surface map presented in Fig. [10] demonstrates that the QWS and the topo- 
logical state undergo an similar hexagonal warping. Both show weaker warping occurring at 
smaller momentum, as described by the model in [73]. While the inner branch of the QWS 
leads to a circular Fermi contour, the outer one appears hexagonal. 

Recently, it has been shown that the dosing of H 2 leads to similar effects as observed for 
CO [42]. The band bending is similar but its origin has been ascribed to a chemical reaction 
occurring at the surface that extracts Se atoms, leaving positively charged vacancies. In the 
case of CO adsorption, such a reaction appears unlikely. CO is a strongly bonded molecule 
(triple-bond), and it can only be dissociated by adsorption on reactive transition metal 
surfaces. The adsorption of CO as a molecular species is consistent with the behaviour of 
the surface electronic structure upon mild annealing of the sample (see Fig. [To] (d)). When 
we heat the sample to 250 K, the band bending and Rashba-splitting is reduced, consistent 
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with a partial desorption of CO. Note that the thermal desorption of CO from Bi 2 Se3 has 
not yet been studied in any detail, and it is not clear if the adsorption can be completely 
reversed by thermal desorption. A recent theoretical study of molecular adsorption of other 
species on Bi 2 Se 3 suggests that this might not be the case due to preferential and strong 
bonding to defect sites on the surface [84] . 

We now return to possible alternative explanations of the phenomena discussed here. 
Different interpretations have been proposed, in particular to explain the Rashba-split states 
in the CB region but also for the M-shaped states in the VB region. Wray et al show that 
states of similar appearance can be found in a slab calculation of the electronic structure 
of Bi 2 Se 3 when the energies of the states near the surface are rigidly shifted [39j. Such 
a procedure can be viewed as a crude implementation of a band bending and is thus not 
dissimilar to the picture given here. In fact, it is otherwise very difficult to accurately 
account for the band bending in a supercell-type density functional theory calculation. The 
reason is that the the slab used in such a calculation usually contains a limited number 
of layers, typically between 12 quintuple layers (114 A) [T7j and 25 quintuple layers (240 
A) [85]. Band bending effects that extend to « 300 A into the sample [34], are thus difficult 
to describe accurately with the limited slab size. 

An alternative scenario is discussed in Refs. [1I614I8] . There it is argued that spectroscopic 
changes similar to those reported here could be caused by an increase in the van der Waals 
gap(s) between the near-surface quintuple layers. This increase of van der Waals gap sizes 
could be caused by the intercalation of adsorbed atoms or molecules. Strongly increased dis- 
tances between quintuple layers (between 10% and 40%) are predicted to lead to two major 
changes in the electronic structure: the appearance of Rashba-split states that are split-off 
the conduction band minimum and M-shaped states in the projected band gap below the 
valence band maximum (see Fig. [jjc)). It is indeed plausible that at least alkali atoms can 
be intercalated into the van der Waals gap in Bi 2 Se 3 , as this is a well-known phenomenon 
for bulk crystals [86], but for molecules such as CO intercalation appears unlikely. On the 
other hand, this proposed mechanism fails to explain a key experimental observation; a cal- 
culation with an increased interlayer spacing predicts a single M-shaped state in the valence 
band region whereas several such states are observed experimentally. This shortcoming can 
probably be improved by allowing for several van der Waals interlayer spacings to relax [48] . 
A more important issue than the number the theoretically predicted M-shaped states is that 
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their energy position of does not agree with the experimental result. In the calculation, the 
M-shaped states are found in the projected band gap below the top of the valence band 
at f . These states are thus genuine new surface states. In the experiment, on the other 
hand, the M-shaped states are found within the region of the projected valence band, with 
the possible exception of the state at highest binding energy. This is directly evident from 
comparing the data shown in Figs. [7] and [9| These states, therefore, cannot be surface 
states, but must be QWS derived from the bulk states. Another interesting point is that the 
calculated surface electronic structure is only significantly modified for rather large changes 
in the van der Waals interlayer spacing, on the order of 10 % or more. Experimentally, the 
intercalation-induced changes are often found to be rather small. The intercalation of sub- 
stantial amounts of copper into bulk Bi 2 Se 3 , for example, leads to a change of the van der 
Waals gap of only 3% [87] . Currently, little is known about intercalation-induced changes of 
the near-surface structure of Bi 2 Se 3 ,but such changes could be determined by surface x-ray 
diffraction, at least in principle. 

The band bending picture for the explanation of the M-shaped states in Bi 2 Se 3 has been 
questioned by recent results for the similar compound Bi 2 iiTe 3 [8T] . In Bi 2 Te 3 , similar M- 
shaped and distinct states in the valence band region are found upon exposing the sample 
to air or 7V 2 , but the authors of Ref. [HI] argue that band bending cannot be the origin of 
these structures for two reasons: first, the structure is M-shaped and thus the effective mass 
changes from electron-like in the centre to hole-like on the wings. Thus, it is impossible for 
the entire state to be confined. Second, the band bending for confining the states would 
have to be larger than the total band width, which the authors determine to be around 
300 meV. We believe these arguments are invalid. Concerning the first argument, the shape 
of the state in k\\ is not relevant for the confinement that takes place in the z direction 
perpendicular to the surface. Along k Z) the most advanced GW-calculations show no large 
dispersion of the band, and it is not clear if it is hole-like or electron-like [88]. Concerning 
the second argument, the relevant bandwidth is also that in k Z) not in k\\. Again, GW- 
theory predicts the width of the upper valence band to be very small [88J, and even the 
band structure calculation shown in the supplementary material of Ref. [81 J appears to give 
a large band width only because two bands are involved. Finally, it is not strictly necessary 
to have a band bending that exceeds the total band width in order to confine the valence 
band states. In a situation with a projected band gap below the valence band, the states 
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with the lowest energy can already be confined by a much smaller band bending. A band 
bending larger than the total width is only required if the entire band is to be quantised. 



CONCLUSION 



In conclusion, we have presented ARPES results for clean and adsorbate-covered Bi 2 Se 3 . 
Different bulk dopings were studied. For carbon monoxide or rest gas adsorption, a strong 
downward bending of the near-surface bands was observed. This was found to lead to 
new spectral features in the energy region of both the valence band and the conduction 
band. These changes were explained by a confinement of the states and the formation of 
quantum well states. This interpretation is supported by simple calculations of the band 
bending effects and by photon energy-dependent ARPES measurements that show how the 
band dispersion in the direction perpendicular to the surface is lost and non-dispersing, 
sharp states appear within the energy regions of the projected bulk bands. For strong band 
bending, the conduction band quantum well states are strongly Rashba split. 
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